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In order to detect ultrahigh-energy cosmic rays (UHECR) , atmospheric fluorescence 
light from the trajectory of the extensive air shower may be measured by mirror- 
photosensor systems. In this type of experiment the photon yield from electrons 
exciting air of various densities and temperatures is most fundamental information 
for estimating the primary energy of UHECR. An experiment has been undertaken 
using a 90 Sr /? source to study the pressure dependence of photon yields, and the 
life times of the excited states, for radiation in nitrogen and dry air. The photon 
yield between 300 nm and 406 nm in air excited by 0.85 MeV electrons is 3.73±0.15 
(±14 % systematic) photons per meter at 1000 hPa and 20 °C. The air density and 
temperature dependence is given for application to UHECR observations. 
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1 Introduction 



The fluorescence technique was first successfully used by the Fly's Eye detector 
to explore cosmic rays in the ultra-high energy region [1] . From 1998 the HiRes 
detector, a successor to the Fly's Eye, has been in operation with improved 
resolution in energy, arrival directions and in meaurements of the longitudinal 
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shower development [2]. Fluorescence detectors will be constructed along with 
the surface array at the Pierre Auger Observatory [3] and one of them is now in 
operation in Malargue, Argentina. In these experiments an atmosphere of quite 
low vapour pressure (i.e. dry) is used as a vast scintillator. In contrast, in the 
case of a satellite-based telescope viewing downward (EUSO) [4], fluorescence 
light from extensive air showers is observed mainly over the ocean where air 
is not dry. 

While the fluorescence efficiencies of electrons in the troposphere are funda- 
mentally important, there are only a few measurements which can be applied 
to the UHECR experiments. Bunner summarized in his Ph.D thesis [5] the 
fluorescence radiation from air and showed that photons from the 2nd posi- 
tive (2P) band of the nitrogen molecule and the 1st negative (IN) band of the 
nitrogen ion were significant. Since there are many vibrational states, there 
are many lines between 300 nm and 400 nm. Photon yields listed in Bunner's 
thesis are based on those excited by low energy electrons, deuterons and al- 
pha particles stopped in air (i.e. total absorption conditions) and the errors 
of each experiment are about 30 %. It is important to measure the photon 
yields excited by electrons in various bands under thin target conditions. In 
order to apply them to air shower observations, their density and temperature 
dependencies must be known in each band. 

Kakimoto et al. [6] measured photon yields in thin targets and showed that the 
yields were proportional to the particle energy loss per unit length, dE / dx, 
between 1.4 MeV and 1000 MeV. However, their measurements were limited 
to three main wavelength bands with central wavelengths of 337 nm, 357 nm 
and 391 nm. The results of Kakimoto et al. are about 20-40 % larger than 
those summarized by Bunner at these three wavelengths. They also measured 
the photon yields using the HiRes optical filter which is sensitive in a wide 
band between 300 and 400 nm. Using this wide band measurement, Kakimoto 
et al. estimate the absolute values of the lines not individually measured, and 
they get the values lower than those determined by Bunner. 

Given these previous measurements, it is clear that precise measurements at 
wavelengths other than the three major lines above are required in thin targets, 
since optical filters between 300 nm and 406 nm are used in most cosmic 
ray experiments. We measured the photon yields of electrons in dry air with 
various pressures to apply to the atmosphere in a desert. The photon yields 
were measured through filters of central wavelengths of 314.7, 337.7, 356.3, 
380.9, 391.9 and 400.9 nm and 10 nm bandwidth. Yields in damp air are now 
under measurement and will be reported in a subsequent paper. 

Recently, the HiRes group reported measurements of the energy spectrum of 
UHECRs [7] . There are significant discrepancies in the energy spectrum at the 
highest energies as measured by two experiments, HiRes and AGASA (Akeno 
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Giant Air Shower Array) [8]. Using the values in the present measurement, 
we discuss their effect on the energy determination of UHECR's by the fluo- 
rescence technique. 



2 Experiment 

2.1 Experimental arrangement 



photomultiplier tube 




Fig. 1. Schematic drawing of the chamber (top view). Three PMTs are mounted 
on two sides and the top of the chamber, and they view fluorescence light through 
quartz windows. Optical filters are mounted between the PMTs and the windows. 
Electrons from 9 °Y^ 90 Zr are beamed and detected by a scintillation counter. 

We chose a photon counting and thin target technique to measure the pressure 
dependence of photon yields (the number of photons produced by electrons 
per meter of travel) from nitrogen and air excited by electrons, following the 
method employed by Kakimoto et al. [6]. The cubic chamber used is shown 
schematically in Fig. 1. Three photomultiplier tubes of 2 inch diameter (Hama- 
matsu photon counting H7195PX, named PM1, PM2 and PM3), which were 
selected for low noise, were mounted on two sides and the top of the chamber, 
and view the light through quartz windows. The window material of PMT 
is borosilicate glass and its typical quantum efficiency is drawn by a dashed 
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curve in Fig. 3. The PMT high voltage and gain used in the measurements 
are about 1950V and 6xl0 6 , respectively. A central region (diameter 2.5 cm) 
of each PMT was used for photon counting. Electrons with a maximum en- 
ergy of 2.28 MeV from (5 decay of 90 Sr^ 90 Y^ 90 Zr (3.7 MBq) were beamed 
by a collimator. The energy of each electron was measured by a scintillation 
counter (electron counter). The coincidence between the electron counter and 
any one of the photon counters was used to generate a gate for the ADCs 
and a start signal for a TDC. A beam length of 5 cm was visible to the three 
photon counters and the triggered rate of electrons was about 2.7xl0 3 s" 1 in 
vacuum and 1.5xl0 3 s _1 at 1000 hPa. The electron energy was reduced by 
about 0.17 MeV by a mylar window and the aluminum cover of the source. 

The energy spectrum of the triggered electrons is shown in Fig. 2. As seen 
there, the measured spectrum fits well to the expected decay spectrum of (3 
rays depicted by the solid curve. The threshold and average electron energies 
are 0.3 and 0.85 MeV respectively. The average energies below and above 0.85 
MeV are 0.55 MeV and 1.2 MeV, respectively. The difference in photon yields 
between these two energy regions will also be presented. 

In the present measurements, six narrow band interference filters with central 
wavelengths of 314.7, 337.7, 356.3, 380.9, 391.9 and 400.9 nm were used. The 
bandwidth of each filter was about 10 nm. Hereafter we refer to 316, 337, 358, 
380, 391 and 400 nm, respectively, to represent the main line in each filter 
band. Fig. 3 shows the transmission coefficients of the adopted filters as a 
function of wavelength (provided by the manufacturer, Hi- Technology Inc). 
Vertical lines show the relative photon yields of fluorescence lines in air at 
1000 hPa listed by Bunner [5], where solid and dotted lines correspond to the 
yield from 2P and IN bands, respectively. There remain some lines which are 
not covered in the present experiment, but their contribution to the total yield 
may be less than 5 %. 

In any one run under a particular set of conditions, the number of incident 
electrons at the electron counter was registered. If there was a coincidence in 
150 ns between the electron counter and any one of three photon counters, 
the ADC values of the electron counter and the corresponding photon counter 
were recorded. The time difference between the photon signal and the electron 
signal (which was delayed by 180 ns) was also recorded by the TDC in units 
of 0.5 ns. 

Examples of the ADC and TDC distributions are shown in Fig. 4. The coin- 
cident time of the electron counter and the photon counter is shown as the 
369th channel of the TDC. Since the TDC is started by the photon signal 
and stopped by the electron signal (delayed by 180 ns), the more the channel 
decreases the more the photon signal is delayed. The time resolution of the 
system is measured by using coincident signals of cosmic ray muons between 
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Fig. 2. Energy spectrum of the electrons. Dots represent the observed spectrum and 
the solid curve is that expected from the (3 decay of 90 Y. The threshold and average 
electron energies are 0.3 MeV and 0.85 MeV respectively. 




300 320 340 360 380 400 

Wavelength (nm) 



Fig. 3. Typical example of quantum efficiency of PMT used ( a dashed curve) and 
transmission coefficients of interference filters used in the present experiment. The 
relative intensities of fluorescence lines at 1000 hPa from the nitrogen molecule 
(solid lines) and ion (dotted line) in air from Bunner [5] are also shown neglecting 
their spectral widths. 

PM2 at the top of the chamber and the electron counter, and is found to be 
1.26 ns. Signal can clearly be separated from the background in the TDC dis- 
tribution. Signals above the 369th channel in the lower plot of Fig. 4 are due to 
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Nitrogen 250 hPa 





Fig. 4. An example of ADC and TDC distributions (for 250 hPa, N 2 ). The upper 
plot is a typical single photon pulse height distribution and the lower plot shows 
the time difference between incident electrons and emitted photons. Each channel 
corresponds to 0.25pC and 0.5 ns, respectively. The TDC is started by the photon 
signal and is stopped by the electron signal which is delayed by 180 ns. 

background photon signals prior to the electron signal, which were triggered 
due to the 20 ns pulse width of the discriminator output. 



2.2 Backgrounds and the difference from our past measurements 



There is a background of single photon PMT noise, which is uniformly dis- 
tributed in the TDC distribution as shown in Fig. 4. With a shutter was 
placed in front of each photon counter inside the chamber, some signals were 
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observed even in vacuum runs. Their rates were (0.082 ± 0.06) min -1 when a 
lead brick was not set around the source in Fig. 1 and were reduced to | with 
the lead brick in the present measurement. This background can, according 
to our calculations, be wholly attributed to Bremsstrahlung photons from the 
aluminum foil (0.1 mm thickness) of the source window. There are coincident 
muons measured by the photon counter at the top of the chamber (PM2) and 
the electron counter. This rate is measured to be (0.133 ± 0.008) min -1 . 

The differences of past measurements reported at the 27th ICRC (Ml) [9] and 
the present ones (M2) are summarized as follows; 

(1) The visible beam length was shortened from 10 cm to 5 cm, so that the 
average photon incident angle was reduced from 9.3° to 5.2°. Though we 
corrected the interference filter transmission coefficients for these photon 
incident angles in Ml, the correction has not been applied to M2. The 
differences between both results are within 5%, irrespective in any filter 
band. 

(2) The background from Bremsstrahlung photons around the source region 
in Ml has been reduced by about | with the 5 cm thick lead brick around 
the source in M2. 



2.3 Derivation of photon yields 



The number of signal counts was obtained by subtracting backgrounds deter- 
mined from the background portion of the TDC distribution. This background 
portion can be discriminated from the signal portion as shown in Fig. 4. The 
photon yield per unit length per electron e, is determined as the number of 
signal counts iV divided by the product of the following: the total number of 
electrons /, the length of the fluorescence portion a, the solid angle of the 
PMT Q, the quartz window transmission rj, the filter transmission /, and the 
quantum efficiency QE and the collection efficiency CE of the PMT. 

N 

e ~ Jxaxfix^x/xQExCE ' ^ 



In the following analysis we use a value of / appropriate for the main line in 
each filter band pass. The effect of another line in one filter will be discussed 
later. / is about (4 ~ 5) x 10 s from about 80 hours in each run. Finally, e in 
vacuum run is subtracted from each e determined above and the corrected e 
is determined. An example of typical data sample is shown in Table 1. 
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Table 1 

An example of typical data sample. 



item 


Nitrogen 


Air 


pressure (hPa) 


1000 


1000 


filter (nm) 


337 


337 


observation time (sec) 


284778 


268919 


total number of electrons 


4.317xl0 8 


4.106x10 s 


coincident counts 


34688 


8393 


counts in signal portion 


28004 


3779 


signal counts 


26920 


3148 


photon yield (/m-electron) 


8.38±0.15 


1.094±0.025 


photon yield in vacuum (/m-electron) 


0.03±0.00 


0.028±0.001 


corrected photon yield (/m-electron) 


8.35±0.15 


1.066±0.025 



2.4 Systematic errors 

The systematic errors of the present experiment are summarized in Table 2. 
The main uncertainty is due to the PMT calibrations (QE and CE), which were 
provided by the manufacturer (Hamamatsu Photonics). The CE is defined by 
Hamamatsu as the counts at the anode divided by the number of photo- 
electrons emitted from the photocathode. In this measurement the factory 
used a 25mm diameter area centered on the photocathode, at a PMT gain 
setting of 5xl0 6 and with a discriminator setting at | of the single photon 
peak. These conditions are also used in the present experiment. The number 
of photo-electrons was estimated in DC mode in a separate measurement. 

The filter transmission coefficients are provided by the manufacturer (Hi- 
Technology Inc.) and the error is mainly due to the various incident angles of 
the photons at the filter. We have estimated this value from the results of Ml 
and M2 described in the Section 2.2 (1). The error due to a possible contami- 
nation of some lines which fall in the tail of a filter's transmission can not be 
derived accurately without absolute photon yields of those lines. Referring the 
table summarized by Bunner we assumed these values as square root sum of 
six filter bands, a few % each (a/6 x 2%), and is listed in the fourth line. Other 
systematic uncertainties appeared in Eq.(l) are as follows: I is negligible; a, 
2%; f2, 3%; and rj , 2%. These are added in quadrature and the result is listed 
in the fifth line. 

It should be noted that the temperature has been maintained within ±0.5°C, 
the pressure within ±0.5hPa and the high voltage within ±0.05% in each run. 
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The uncertainty due to the variation of these factors is negligible in estimating 
e, but about a few% in evaluating p' of 391nm line. Therefore total systematic 
error in e is estimated to be 14 %, assuming the errors add in quadrature. The 
statistical error in each run is less than 3%. 
Table 2 

Systematic errors of the present experiment 



No. 


item 


errors 


1. 


Quantum efficiency of PMT 


5 % 


2. 


Collection efficiency of PMT 


10 % 


3. 


Transmission coefficient of filter 


5 % 


4. 


Contamination from lines at the tail of filter transmission 


5 % 


5. 


Other parameters (I, a, f2, rj of Eq.(l)) 


4 % 




Total 


14 % 



3 Results of photon yields 

3. 1 Nitrogen 

As a first step we have measured e in nitrogen, so that we may compare results 
with previous experiments. The pressure dependencies of e in six wave bands 
are indicated by solid circles in Fig. 5. In each figure, the main line in each 
filter band is indicated. The results of Kakimoto et al. (cited from Ueno [11]) 
are also plotted with open squares in the 337 nm and 358 nm bands. The 
present results at 337 nm and 358 nm are in good agreement with those by 
Kakimoto et al., even though the PMTs, filters and gas chambers are different 
in the two experiments. 

The solid lines in the figure are from Eq.(10), which will be discussed in Section 
4.1. A dotted line in the plot for the 391 nm band is the best fit to the equation 
with the five highest pressure points excluded. The experimental points in 
some filters cannot be fitted with a constant p' throughout the pressure range, 
especially in the 391 nm band. This may be due to the superposition of two 
lines in one filter band, as will be discussed in Section 4.3. 

3.2 Air 

We used a mixture of 78.8 % nitrogen gas and 21.2 % oxygen gas, without 
Argon and the other small contributions to air, such as C0 2 , Ne, CO, water 
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Fig. 5. The pressure dependence of e in N2 at 17 °C excited by electrons with an 
average energy of 0.85 MeV in six wave bands. In each figure, the main contributing 
line in each wave band is listed. Solid circles are the present results. Open squares 
are from Kakimoto et al. cited from Ueno [11]. Solid curves are the best fits of 
Eq.(10) with p' as shown, discussed below in Section 4.1. The dotted line for the 
391 nm band shows the best fit when the five highest points are excluded. 



vapour etc. This air-like mixture is expected to give similar photon yields to 
dry air, since the change of photon yields may be negligibly small even if Argon 
of 1 % is included [5]. In Kakimoto et al, a mixture of nitrogen and oxygen 
gas was also used as dry air. 

The pressure dependence of e in dry air is shown in Fig. 6 by filled circles 
in six wave bands. The results of Kakimoto et al. [6] are also shown by open 
squares in the 337 nm, 358 nm and 391 nm bands. Their results are larger 
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Pressure (hPa) Pressure (hPa) 

Fig. 6. The pressure dependence of e in air at 20 °C. The data of Kakimoto et al. 
in dry air at 15 °C with 1.4 MeV electrons is plotted by open squares. Solid lines 
show the best fit of Eq.(10) with the value of p' as shown, as discussed in Section 
4.1. The dotted line in the plot of the 391 nm dependence is the best fit excluding 
the four highest points. 



than the present ones by about 0.1 photons/m irrespective of the pressure, and 
in all three bands. This is understood as follows: there were backgrounds of 
photons from Bremsstrahlung at the foil of the source window in Kakimoto et 
al., while those photons were reduced by a lead brick in this experiment. The 
effects of this background were negligibly small in the nitrogen measurements 
as shown in Fig. 5. The solid lines come from Eq.(10), which will be discussed 
in Section 4.2. 
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4 Analysis 



4-1 Fluorescence decay time and efficiency in N 2 gas 



The observed decay time r of the fluorescence is related to the following three 
terms: the mean lifetime of the excited state for decay to any lower state r r , 
the lifetime of collisional de-excitation r c and that of internal quenching 7$. 
This can be expressed as follows [5], 

111111 

- = - + - + - = - + -, (2) 



and 



V^MkTl 

r c = — : , 3 

4o- n „ p 



where M is the N 2 molecular mass, k is Boltzmann's constant, T is the gas tem- 
perature, o~ nn is the cross-section for nitrogen-nitrogen collisional de-excitation 
and p is the gas pressure. We define the reference pressure p' [5] , which cor- 
responds to the pressure where r is equal to r c , as 

. VnMkT 1 

P = ^0—7 ■ (4) 

^ u nn 'o 



Then the reciprocal of the lifetime ^ is related to p by 

-= \^-]p+- = -(ap+l) , (5) 



where a = -7. 

v 

The fluorescence efficiency of the i-th band at pressure p, is defined as 

the power radiated by the gas in the i-th band per unit power deposited in 
the gas by an electron. This can be written as the ratio of de-excitation 

by radiation to the total, and is expressed as a function of p as follows [5] : 

Since the energy available to produce photons in the i-th band is the electron 
energy loss per unit length times $i(p), the photon yield per unit length per 
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electron, for i-th band under atmospheric pressure p is written as follows 
[10]: 



and 



Hp) = 



(8) 



where p is the air density (kg m 3 ), — — is the energy loss in (eV kg 1 m 2 ), 

and hvi is the photon energy of the i-th band (eV). $° (=^) corresponds 
to the fluorescence efficiency for the i-th band in the absence of collisional 
quenching. 

By using the equation of state of a gas, p = pR^ 2 T where R^ 2 is the specific 
gas constant in N 2 (296.9 m 2 s _2 K _1 ), Eq.(7) is rewritten as a function of p at 
a constant temperature T in Kelvin: 




(9) 
(10) 



where x — - and a, — \. This relation can be written in a more general form 

p Pi 

as a function of the gas density and the temperature as [6] 



where 



^ = (S)fi and B = = ^ToRy, 

hvi p\ y/nkM 



The pressure dependencies of - for various wave bands in N 2 are shown in 
Fig. 7. Each r and its standard deviation in the figure was calculated with 
a weighted least square by fitting the counts(J) in the decay portion of the 
TDC distribution (see Fig. 4) to the equation In J = In J — K Due to the time 



13 




Fig. 7. The pressure dependencies of ^ for the excited states of N2 for 6 wave bands. 

resolution of the present experiment (1.26nsec) and the limited statistics, the 
error in ^ of each point is quite large compared with that of e. 

In order to fit the experimental points of Fig. 5 and 7 to Eq.(lO) and (5) 
respectively, the least square (LS) method is utilized. Let us consider first the 
fit to the p — e diagram (Fig. 5), since the data points are more numerous than 
those in p — ^ diagram (Fig. 7) and the large difference of errors between both 
diagrams as mentioned above. For each i-th band the first x\ is defined by 

^5>fe-^) 2 , (13) 

3 J 



where the weights Uj come from the errors in €j at Pj(— ^7). By minimizing 

xj, we find the first values of and Then after linearizing Eq.(lO) 

using and Newton's method [13] is applied to estimate the final 
and after the z/-th iteration. The statistical errors of and are 
evaluated using the propagation law of errors using those of ej. For a given 
a,("\ it is easy to find r and its error by the LS method because Eq.(5) is 
linear with respect to (ap + 1). The obtained results are shown by solid curves 
in Fig. 5 and by lines in Fig. 7. The parameters are listed in Table 3. The 
values of a nn derived from tqp' and $° from C are also listed in the table. In 
the case of the 391 nm band, we find that the experimental points do not fit 
well to one value of p' throughout the pressure range but a dashed line in Fig. 
5 shows Eq.(lO) for a limited range of pressures. 
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Table 3 

t , p', cr n n, C and <£>° in N 2 . 



wavelength 


To 


p> 


"ran 


C 


$° 


nm 


ns 


hPa 


xlQ -20 m 2 


xl0~ 2 /(hPa-m) 


xl0~ 3 


316 


20.5± 


2.2 


88.3 ± 7.5 


3.36± 0.46 


2.50±0.14 


0.509±0.028 


337 


26.8± 


1.4 


151.8 ± 4.8 


1.50± 0.09 


6.27±0.13 


1.20 ±0.03 


358 


30.2± 


0.6 


118.0 ± 4.0 


1.71± 0.07 


5.71±0.13 


1.03 ±0.02 


380 


30.9± 


0.6 


115.9 ±11.1 


1.70± 0.17 


2.26±0.10 


0.382±0.017 


391 


23.4± 


0.6 


4.42± 0.28 


59.0 ± 4.1 


11.5 ±0.5 


1.90 ±0.09 


400 


23.8± 


1.8 


94.2 ± 7.9 


2.72± 0.31 


0.72±0.04 


0.116±0.007 



Table 4 

Comparison of $j(p) in N 2 at p=800 hPa (600 mmHg). D. and O. refers to Davidson 
and O'Neil [14k 





This experiment 


Ueno 


D. & 0. 


beam 


electron 


electron 


electron 


energy 


0.85 MeV 


1.4 MeV 


50 keV 


wavelength 


xlO" 4 


316 nm 


0.505±0.051 






337 nm 


1.90 ±0.07 


1.87±0.01 


5.20 


358 nm 


1.32 ±0.05 


1.35±0.01 


3.70 


380 nm 


0.483±0.051 




1.40 


391 nm 


0.104±0.007 




0.102 


400 nm 


0.122±0.013 




0.200 



Our values of in N 2 at 800 hPa are listed in Table 4, together with those of 
other experiments. The electron energy loss 4^ is calculated using the equation 
in Sternheimer et al. [15] with a density correction, and is 0.1687 MeV/kg-m~ 2 
and 0.1677 MeV/kg-m~ 2 at 0.85 MeV, for N 2 and air respectively. The values 
from the experiments by Kakimoto et al. are cited from Ueno [11] and agree 
well with the present ones for the 337 nm and 358 nm bands. However, the 
values by Davidson & O'Neil [14] are about a factor of 2.8 larger than the 
present results except at 391 nm and 400 nm. 

The present result for $° at 391 nm is (1.90±0.09)x 10~ 3 which does not 
agree with the result (6. Oil. 9) x 10 -3 by M.N.Hirsh et al. [12], measured 
using electrons of 1.46 MeV under a gas pressure of less than 10 hPa. This 
will be discussed in Section 4.3 in connection with the two lines contained 
within this filter band. 
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Fig. 8. The pressure dependence of the reciprocal of the fluorescence decay time (^) 
in air. 

4-2 Fluorescence decay time and efficiency in dry air 



In the case of air, the reference pressure p' of Eq.(4) may be written as [5]: 



4r 



p' ^irM n kT 



fnO~nn f o0~r, 



2M n 



D 

7t ' 



(14) 



where M n and M Q are the masses of nitrogen and oxygen molecules, respec- 
tively. a nn and a no are the cross-sections for collisional loss by excited nitrogen 
molecules without radiation, through collisions with other nitrogen or oxygen 
molecules, respectively. The fraction of nitrogen f n is 0.79 and that of oxygen 
fo is 0.21. 

The pressure dependence of ^ in dry air is shown in Fig. 8 for six wave bands. 
The fitting of parameters from Eq.(10) and (5) to the data points in Fig. 
6 and 8 is done using the same method as described for the N2 case. Three 
parameters, a (or p'), C and r G and their errors are thus evaluated. The results 
are shown by solid curves in Fig. 6 and by lines in Fig. 8 and the parameters 
are listed in Table 5. Since p' is determined from the pressure dependence of 
e, the line of 391nm in Fig. 8 seems not to be best fitted to the experimental 
points. The cross-section a no is derived from r p', where p' is given by Eq.(14) 
with <7 nn determined previously. $° and C are also included in the table. In 
the case of 391 nm, the dashed curve in Fig. 6 represents parameters in Table 
5. 
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Table 5 

r , p', a no , C and <&° in air 



wavelength 


To 


v' 


""no 


C 


3>° 


nm 


ns 


hPa 


xl0 -19 m 2 


xl0" 2 /( nPa - m ) 


xl0~ 3 


316 


22.0±0.5 


24.0 ±2.0 


4.30±0.70 


2.34±0.14 


0.463±0.029 


337 


32.4±1.8 


19.2 ±0.7 


4.06±0.29 


5.43±0.15 


1.01 ±0.03 


358 


34.2±0.4 


19.3 ±1.4 


3.83±0.31 


4.86±0.27 


0.850±0.047 


380 


33.9±1.6 


21.2 ±2.6 


3.45±0.46 


1.76±0.18 


0.289±0.029 


391 


24.2±3.2 


4.61±0.34 


3.70±0.62 


7.04±0.39 


1.13 ±0.06 


400 


26.9±1.9 


19.3 ±1.8 


4.65±0.76 


0.86±0.06 


0.134±0.009 



Table 6 

Comparison of $j(p) in air at p=800 hPa 





This experiment 


Kakimoto 


D. & O. 


Hartman 


Bunner 


beam 


electron 


electron 


electron 


electron 


a 


energy 


0.85 MeV 


1.4 MeV 


50 keV 




4 MeV 


nm 


xl0~ 6 


316 


13.1 ±1.4 






3.0 


5.3 


337 


23.7 ±1.1 


21 


21.0±3.2 


20.0 


8.5 


358 


20.2 ±1.8 


22 


15.0±2.3 


16.0 


5.8 


380 


7.5 ±1.2 




5.2±0.8 




4.1 


391 


6.48±0.60 


8.4 


7.0±1.1 


5.9 


4.8 


400 


3.18±0.36 




1.8±0.3 


1.7 


1.5 



The yield in air is expressed by Eq.(ll), after replacing R N2 by i? air =287.1 
m 2 s _2 K _1 and by using p\ given by Eq.(14). The values of in air at 800 
hPa are listed in Table 6 together with results of other experiments. Those of 
Hartman and Bunner are calculated by using <3>° and p' listed in Table 2.6 of 
Bunner [5]. These values of at p=800 hPa are in good agreement, except 
those of Bunner which were measured using alpha beams. 

In order to see the energy dependence of in air, values of e« are derived 
in two energy regions below and above the average energy (0.85 MeV). In 
Fig. 9, the energy dependence of q in air at 1000 hPa is shown for two wave 
bands together with measurements by Kakimoto et al. cited from the thesis 
of Ueno [11]. The of an electron in air is calculated from a formula given 
in Sternheimer et al. [15] and the solid line is normalized to the 1.2 MeV yield 
point of the present measurements. A dashed line is normalized to the 1.4 
MeV point of Ueno in each figure. The solid lines agree with the higher energy 
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Fig. 9. Energy dependence of e in air at 1000 hPa. The dE/ dx curve is shown by 
a solid line, and is normalized to the yield at 1.2 MeV (present measurement). 

measurements (between 300 MeV and 1000 MeV) better than the dashed 
lines. This may be because the photon backgrounds were not excluded from 
the 1.4 MeV points of Ueno, as described in Section 3.2. Our measurements 
are consistent with e being proportional to 4p. 



4-3 Two components analysis in one filter band 



There are some discrepancies between the experimental points and the model 
pressure dependencies of ^ and q as shown in Eqs.(5) and (10), especially in 
the 391 nm filter band. In the following, we try to fit the observed pressure 
dependencies of — and e, with a superposition of two lines in one filter band. 
In this case the observed photon yield e G bs(p) is the sum of the photon yields 
of the main line e\{p) and the sub-line e 2 (p), and is written as follows by 
extending Eq.(10): 



C C 

e Q bs(p) = ei(p) + e 2 {p) = — ^— + — — — , (15) 

x + ai x + a 2 

where x — ^, a± — ^- and a 2 — C\ and p[ are the parameters of the main 
line in the filter band, and C 2 and p' 2 are the parameters of another line. It 
should be noted that e Q b s is determined assuming the filter transmission of the 
main line in the filter, and hence e 2 must be corrected with using the filter 
transmission of that line. 
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The reciprocal of the observed life time - is approximately expressed by 

the weighted mean of ^-^y and ^^y- The weights are expressed in terms of 
the relative photon intensities of two lines. Therefore, 



1 _ ( \ 1 , ( € i(p) \ 1 ( 16 ) 



^obs (p) V e i (?) + e 2 (P) J r i (p) V e i (p) + e 2 (p) / r 2 (p) 
where 



1 1 / 

-(aip + 1) , 



n(p) r i 

J_ = J_( a2 p + l) . (17) 

T 2 (p) T o2 



In this case, we have to determine a set of four parameters a 1; Ci, a 2 and 
C 2 in Eq.(15). For many such sets, the value of x 2 is computed and the set 
having the minimum \ 2 is taken as the initial set. Newton's method is also 
applied for this set, which leads to the final set after the v-th iteration. Values 
of a±(— 4-), Ci, a 2 (= -V) and C 2 determined in this way are listed in the upper 

half of Tables 7 and 8. The values of minimum x 2 are also depicted as x 2 i n 
the tables. 

The errors of these four parameters (as listed in the tables) are estimated in 
the following way. Since the probability density function (pdf) of A% 2 , i.e. the 
deviation from Xmm> follows the pdf of x 2 with four degrees of freedom, we take 
AXcriticai = 4-72 where the coverage probability is 0.683 corresponding to the 
1-sigma [18]. Many sets of (p'[, C[, p 2 , C 2 ) near to (p[, Ci, p' 2 , C 2 ) are selected 
such that the value of x 2 almost equals Xmin + ^Xcriticai- The maximum and 
minimum values of {p'{ — p[) lead to the error bounds of p[. Errors of other 
three parameters are estimated in the same way. 

From Eq.(15), (16) and (17), - is expressed as 



P xv , (18) 



T"obs(p) ei(p) + e 2 (p) 
where 



v = ^ + ^ . (19) 

Tol T o2 

That is, t i and r o2 can't be evaluated independently without any assumptions 
and only the value v is determined by the LS method. It should be noted that 
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v is almost insensitive to p' and C, since these parameters are determined to 
fit e G 6s- 
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Fig. 10. Two line fitting for the 391 nm bands in N2. The upper figure is the best 
fit combination of two lines. The lower figure is the best fit for a fixed p[. The 
contribution from the main line is shown by a dashed curve, and the secondary line 
by a dotted curve. The solid curves are the sum of two lines in the left-hand figures 
and the solid line is Eq.(18) with v listed in the right-hand figures. 

Table 7 

Parameters of two line fitting in the 391 nm band of N2. 



Wave length 


p> 


V 


C 


<J>° 


x 2 


nm 


hPa 


10 _3 /(hPa-m-ns) 


10" 2 /(hPa-m) 


10~ 3 




391 


3.97±0.37 


4.79±0.36 


11.8 ±0.7 


1.94 ±0.11 


32.7 


394 


qqc: +242 




0.071±0.028 


0.012±0.005 




391 


1.40 


4.79±0.36 


12.0 ±2.0 


1.95 ±0.33 


158. 


394 


12.6 til 




3.38 ±0.26 


0.55 ±0.04 





The best fit result for a superposition of two lines in the 391 nm band of N 2 
is shown in the upper figure of Fig. 10, where the contribution from the main 
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line is shown as a dashed curve and that from the second line is shown as a 
dotted curve. The solid curve in the left-hand figure is the sum of the dashed 
and dotted lines. In the right-hand figure, a solid line is drawn with the value 
v. The parameters of each line are listed in Table 7. 

It is clear that the fitting is much better than the single line fitting in Fig. 5. 
However, the value of <3>° at 391 nm is much smaller than the measurement 
of Hirsh et al.[12], which was determined using 1.45 MeV electrons in the 
pressure range between and 10 hPa. If we fit their pressure dependence of 
photon yields with Eq.(10) and with fixed p[=lA0 hPa, the best fit case is 
shown in the lower half of Fig. 10 and the parameters are listed in the lower 
half of Table 7. The difference of $° at 391 nm is within experimental errors, 
but the difference of $° from Hirsh et al. is still beyond the measurement errors 
of both experiments. To resolve the issue, we need to measure the yields at 
391 nm with a narrower bandwidth filter at low pressure. 





Fig. 11. Two line fitting of the 391 nm band in air. Other details are the same as 
in Fig. 10. 

In the case of the 391 nm band in air, two cases are also shown in Fig. 11. As 
before, the first shows the best fit combination of two lines, allowing freedom 
in the pressure dependence of each line. The second set of plots shows the case 



21 



Table 8 



Parameters of two line fitting in the 391 nm band in air. 



Wave length 


p> 


V 


C 


3>° 


x 2 


nm 


hPa 


10~ 3 /(n pa - m - ns ) 


10- 2 /( hPa - m ) 


1(T 3 




391 


2.45±0.85 


2.89±0.26 


8.34±1.44 


1.33±0.23 


26.8 


394 


21.8 ±8.5 




0.77±0.33 


0.12±0.05 




391 


1.27 


2.90±0.26 


9.98±0.14 


1.60±0.02 


33.0 


394 


13.4 ±2.3 




1.80±0.35 


0.29±0.06 





where we use a fixed value of p[ as determined by Hirsh et al. [12]. 



In the other filter bands we have also determined the best fit p' and r values 
for two lines in a filter band. However, the best fit values of p' are quite different 
in some cases for the 337 nm, 356 nm, 380 nm and 400 nm bands, where lines 
of the 2P(0,f") and 2P(l,t>") transitions are superimposed. We have decided 
not to list the values in this article and we continue the experiments with a 
narrower band width and improvements of the analysis method. 



5 Discussion 

5. 1 Summary of the results 

In Table 9, we summarize the present results for electrons of average energies 
of 0.85 MeV at 1000 hPa. Fitted values of e in the 4th column correspond to 
photons in each filter band, assuming only one line in each band, and assuming 
a filter transmission at the wavelength of the main line. The exception is the 
391 nm filter, where e is separated into two lines. The values in parentheses 
with upper suffix a in column 4 are not measured in this experiment but are 
estimated as follows. The average ratio of our measured values in other wave 
bands to the corresponding values listed in Bunner [5] is determined, and this 
ratio is multiplied by the yields in Bunner. In the case of the filters with 
central values of 380.9 nm and 400.9 nm, the transmission coefficients for the 
375.6 nm and 405.9 nm lines are low compared to the values used and the 
photon yields corrected for transmission coefficients of these lines are listed in 
parentheses and denoted by an upper suffix b . The sum of the measured and 
corrected values are listed at the bottom of column 4. 

Values of $° calculated from Eq.(10) are listed in the 5th column. In order to 
apply the present measurements to air shower fluorescence experiments, the 
fluorescence yields per meter per electron must be known as a function of air 



22 



Table 9 

Summary of the present experiments. Photon yields (e) per meter per electron of 
average energy 0.85 MeV are values in dry air at 1000 hPa and 20 °C. Values in 
parentheses with an upper suffix a have not been measured in this experiment and 
are estimated as described in the text. Values in parentheses with an upper suffix 



are corrected values taking into account the transmission coefficient (see text for 
detail). 



i 


A 


transition 
state 


e 




A 


B 


nm 


m _1 


xl0~ 3 


m 2 kg _1 


m 3 kg _1 K~5 


1 


311.7 


2P(3,2) 


0.552±0.033 


0.463±0.029 


19.6 ±1.2 


2.04±0.17 


2 


313.6 


2P(2,1) 


3 


315.9 


2P(1,0) 


4 


328.5 


2P(3,3) 


(0.041±0.012) a 








5 


330.9 


2P(2,2) 


6 


333.9 


2P(1,1) 


1.024±0.028 


1.006±0.027 


45.4 ±1.2 


2.56±0.09 


7 


337.1 


2P(0,0) 


8 


346.9 


2P(3,4) 


(0.030±0.009) a 








9 


350.0 


2P(2,3) 


10 


353.7 


2P(1,2) 


0.926±0.051 


0.850±0.047 


40.7 ±2.3 


2.54±0.18 


11 


357.7 


2P(0,1) 


12 


367.2 


2P(3,5) 


(0.054±0.016) a 








13 


371.1 


2P(2,4) 


14 


375.6 


2P(1,3) 


(0.493±0.074) fe 
0.367±0.038 


(0.388±0.058) 6 
0.289±0.029 


(19.8 ±3.0 ) b 
14.7 ±1.5 


2.31±0.25 


15 


380.5 


2P(0,2) 


17 


391.4 


1N(0,0) 


0.204±0.079 


1.33 ±0.23 


69.8 ±12. 


20.1 ±6.9 


16 


389.4 


2P(3,6) 


0.163±0.071 


0.121±0.052 


6.4 ±2.8 


2.25±0.09 


18 


394.3 


2P(2,5) 


19 


399.8 


2P(1,4) 


0.163±0.011 
(0.247±0.043) 6 


0.134±0.009 
(0.205±0.035) 6 


7.2 ±0.5 
(10.9 ±1.9 ) b 


2.54±0.24 


20 


405.9 


2P(0,3) 


Sum(measured) 


3.40 ±0.13 








Sum(corrected) 


3.73 ±0.15 









density and the temperature for each wavelength. A^ and Bi from Eq.(ll) are 
listed in the 6th and 7th column. 
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The total photon yield between 300 nm and 406 nm given by Bunner is about 
14 % smaller than the present corrected results. His value was obtained as 
weighted averages of three measurements by Davidson and O'Neil (50 keV 
electrons) [14], Hartman (electrons, cited from [5]) and Bunner (alpha par- 
ticles) [5], and the accuracy of each of the three experiments was not better 
than ±30 % [5] . Therefore the disagreement may be within their measurement 
errors. 



5.2 Effect of the present results on the energy estimation of UHECR 




Fig. 12. Photon yield between 300 and 406 nm for a 0.85 MeV electron as a function 
of altitude. The US Standard Atmosphere 1976 [19] is used for the atmospheric 
model. The thin solid line is the altitude dependence of the 391 nm line and the 
thin dashed lines are those of other lines. The heavy dashed line is the sum of all 
measured lines. A heavy solid line represents the corrected values for the unmeasured 
contributions. The dotted heavy line is the photon yield between 300 nm and 400 
nm from Kakimoto et al., calculated for a 0.85 MeV electron using their equation. 

In Fig. 12, the photon yield between 300 nm and 406 nm for a 0.85 MeV 
electron is shown as a function of altitude. The US Standard Atmosphere 
1976 [19] is assumed. The thin solid line is the altitude dependence of the 391 
nm line and the thin dashed lines are those of the other lines listed in Table 9. 
The heavy dashed line is the sum of all measured lines. It should be noted that 
the altitude dependence is almost the same for all 2P lines. The yield from the 
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391 nm line is, however, almost independent of altitude, but its proportion 
increases with altitude. The heavy solid line is the corrected one described 
in the previous section. As done by Kakimoto et al., the total photon yields 
between 300 nm and 406 nm can be approximated as a superposition of two 
sets of terms, for the 2P lines and the IN line, as follows: 



dE, , ^ A 

. .. ill! V , ^ I 

dE 
dx 



dx 



xp( , + — ' 7= ) , (20) 



)o.85MeV \ 1 + pBiVT 1 + pB 2 \/T 



where p is the air density in kg m 3 and T is in Kelvin (K). Our values are A\ = 
147.4 ±4.3 n^kg- 1 , A 2 = 69.8 ±12.2 n^kg" 1 , B l = 2.40 ±0.18 nr'kg^K- 1 / 2 , 
B 2 = 20.1 ± 6.9 m 3 kg _1 K -1 / 2 . Here the energy loss is normalized at 0.85 
MeV. Though the parameters A 1: A 2 , B x and B 2 from Kakimoto et al. [6] are 
different from the present results, their altitude dependence is quite similar to 
the present experiment, as shown by the heavy dotted line in the figure. But 
the absolute values of the yield are somewhat different. 




Fig. 13. Comparison of e from the present experiment (solid line) with those from 
Bunner (dotted line) and those used in the HiRes analysis (dashed line). 

Fig. 13 shows a comparison of e from the present results between 300 nm 
and 406 nm, with those from Bunner and those values used by the HiRes 
experiment [17]. For an electron in air, (-^f)o.85MeV=0.1677 MeV/kg-m~ 2 and 
(-ff)i.4MeV— 0.1659 MeV/kg-m~ 2 corresponding to ours and the HiRes, re- 
spectively. Therefore we may directly compare the present values of yields 
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and those used by them. Though the values used in the HiRes experiment are 
larger than the present ones in the three main bands (337 nm, 356 nm, 391 
nm), the HiRes total number of photons between 300 and 406 nm is about 13 
% lower than ours as seen in Fig. 13. Though the systematic uncertainty of 
the present experiment is 14%, the main uncertainties are due to QE and CE 
and these measurement procedure by the factory is common to the present 
and Kakimoto et al.. Since the yields used in HiRes code is based on Kakimoto 
et al. [17], the primary energy of UHECR may be estimated larger when we 
use the HiRes photon yields than the present results. In order to know how 
much the energy is largely estimated, it is necessary to take into account the 
wavelength dependence of the attenuation of photons in the air (Rayleigh and 
Mie scattering), the optical filter used, the quantum efficiency of the PMT's 
and other factors which depend on the wavelength. 

Considering only Rayleigh scattering and the optical filter transmission coeffi- 
cient for the HiRes experiment, we have estimated a change in energy of HiRes 
events if we use the present results instead of their values. The HiRes result 
may be largely estimated by about 10 % for showers of a distance around 10 
km and this value increases as the distance increases. We must evaluate the 
factor in more detail, by taking into account the seasonal and daily variation 
of the density, temperature and the Mie scattering at the experimental site, 
and the reconstructed geometry of each shower track. 



6 Conclusion 

Photon yields have been measured in six wave bands as a function of pressure, 
for nitrogen and dry air excited by electrons of an average energy of 0.85 MeV. 
The pressure dependencies of fluorescence decay time were also measured. The 
results are summarized as follows: 

(1) Nitrogen results: The fluorescence efficiencies at 800 hPa are in good 
agreement with the results by Kakimoto et al. [6,11] at 337 nm and 
358 nm, but are smaller by a factor of 2.8 compared with the values of 
Davidson and O'Neil [14] at 337, 358 and 380 nm. On the other hand, at 
391 nm the present result agrees with theirs. 

(2) Air results : The fluorescence efficiencies at 800 hPa are in good agreement 
with the previous experiments [6,11,14,16], if we take into account the 
backgrounds which were not subtracted in Kakimoto et al. However, $° 
for 391 nm is smaller, by a factor 3.5, than the measurement of Hirsh et 
al. [12], provided we use the best fit value of p' from our experiment. 

(3) The photon yield between 300 nm and 406 nm at 1000 hPa and 20 °C is 
3.73±0.15 per meter for an electron of 0.85 MeV. The systematic error 
is 14 %, with the main contribution relating to the collection efficiency 
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of the PMT. The photon yield is proportional to 4-, and its density and 
temperature dependence is given by Eq.(20). 
(4) The photon yield between 300 nm and 406 nm used by the HiRes exper- 
iment is about 13 % smaller than determined by the present experiment. 
If we take into account the wavelength dependence of Rayleigh scattering 
and the transmission coefficient of the HiRes filter, the primary energies 
of UHECRs based on their photon yield is estimated larger than those 
based on the present experiment. The numerical factor depends on the 
distance between the shower trajectory and the experimental site. 
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